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The present study shows the existence of intrinsic coastal air-sea coupled phenomenon in the coastal ocean 
off Baja California and California in boreal summer for the first time. It contributes significantly to the 
interannual sea surface temperature (SST) anomalies there. An initial decrease/increase in the equatorward 
alongshore surface winds weakens/strengthens the coastal upwelling and raises/lowers the coastal SSTs 
through oceanic mixed-layer processes. The resultant coastal warming/cooling, in turn, heats/cools the 
overlying atmosphere anomalously, decreases/increases the atmospheric pressure in the lower troposphere, 
generates an anomalous cross-shore pressure gradient, and thus reinforces or maintains the alongshore 
surface wind anomalies. The regional air-sea coupled phenomenon seems to be analogous to the well-known 
El Nino/Southern Oscillation (ENSO) in the tropical Pacific but with much smaller time and space scales, 
and may be referred to as California Nino/Nina in its intrinsic sense. 
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Climatological surface winds blow southward all year round along the western coast of Baja California and 
California and drive the surface water offshore owing to the Ekman transport' originated in the Coriolis 
force that deflects moving objects to the right of their paths in the Northern Hemisphere. To compensate 
for the surface water, cold and nutrient-rich subsurface water upwells and increases biological productivity in the 
coastal surface layer, which supports abundant marine species at different trophic levels. In some years, however, 
the upwelling is reduced (enhanced) and thus the coastal ocean becomes warmer (colder) than usual 2 " 5 . For 
example, the El Nino years of 1957/58, 1976/77, 1982/83 and 1997/98 witnessed pronounced warming in the 
northeast Pacific; the anomalous warming lasted for over one season with the maximum monthly SST anomalies 
up to 2-6°C 2 " 5 . The anomalous warming caused general decline in biomass of phytoplankton and zooplankton 
and intrusion of warm-water species into higher-than-usual latitudes of habitats 6 " 9 . Extreme climate was also 
observed in the nearby continent; increased winter storms hit California in 1982/83 and resulted in one of the 
wettest winters in record. In contrast, the coastal ocean was extremely cold in the 1 999 La Nina year and California 
experienced cooler-than-normal summer 10 . Researches on the frequent co-occurrences of the tropical ENSO and 
the coastal warming/cooling have shown that ENSO can influence the coastal ocean via atmospheric teleconnec- 
tions that modify strength of the alongshore surface winds 4,11,12 and/or via coastally trapped Kelvin waves that 
propagate poleward from the eastern tropical Pacific and influence the coastal thermocline depth, upwelling and 
SSTs 13 . However, not all of the coastal warming/cooling events are associated with ENSO. The basin-wide 
atmospheric circulation anomalies in mid-latitudes related to the Pacific Decadal Oscillation and the Aleutian 
Low are also regarded as contributors to the coastal SST variations 12,14 ' 15 . 

In general, subtropical oceans are passively forced by the overlying atmosphere and have weak influences on 
the atmosphere, compared to the atmospheric internal variability, in seasonal to interannual time scales 16 . 
However, the SST anomalies in the coastal upwelling ocean may affect the land-sea thermal contrast and 
cross-shore pressure gradient in the lower troposphere, thus influencing the strength of the alongshore surface 
winds. This hypothesis associated with the feedback between the coastal upwelling and alongshore surface winds 
is adopted to partly explain the seasonal intensification of subtropical highs in summer in the Northern 
Hemisphere 17 . Recently, the coastal ocean-atmosphere coupled process is named the coastal Bjerknes feedback 
and applied to explain the regionally intensified SST variability in the coastal ocean off western Australia 18 . Here, 
we demonstrate that the same coastal Bjerknes feedback contributes significantly to the coastal SST variations off 
Baja California and California in boreal summer and identify the existence of a regional coupled mode as the 
California Nino/Nina in relation to the pioneering work of Simpson (1983) among others that describe the 
warming events in the coastal ocean off western North America as California El Nino 319 . The important difference 
from previous works, however, is that we demonstrate the existence of an intrinsic regional coupled mode that is 
independent of ENSO events. 

In the present study, we adopt Monthly National Oceanic and Atmospheric Administration (NOAA) 
Optimum Interpolation Sea Surface Temperature (OISST) version 2 20 , Pacific Fisheries Environmental 
Laboratory (PEFL) derived upwelling indices 21 , and National Centers for Environmental Prediction (NCEP)/ 
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National Center for Atmospheric Research (NCAR) reanalysis l 22 for 
the period of January 1982 to December 2011. Monthly averaged 
anomalies in SSTs, upwelling indices and atmospheric variables are 
calculated by removing the mean seasonal cycles. As shown in Fig. la, 
the coastal SSTs off Baja California show substantial interannual 
variations. The spatial distribution of monthly standard deviation 
of SST is similar to the first mode of empirical orthogonal function 
(EOF) analysis that explains —47% of total variances (Fig. lb). The 
California Nino/Nina index is defined as the SST anomalies averaged 
over 110°W-120°W and 20°N-30°N (the enclosed coastal ocean in 
Fig. la) where they show the largest interannual variability. In the 
past three decades, the California Nino/Nina index (Fig. lc) is 
significantly correlated with the time series of the first EOF mode 
with the correlation coefficient up to 0.9, significant at a 99.99% 



confidence level by the two-tailed r test. Hence, both time series 
represent almost the same interannual SST anomalies off Baja 
California and California. 

Seasonally stratified standard deviations of the California Nino/ 
Nina indices show strong interannual variations from boreal sum- 
mer to winter. As mentioned above, the California Nino/Nina is 
strongly influenced by the tropical ENSO. However, the influences 
are much stronger in boreal winter than summer; the correlation 
coefficient between January-March (July- September, JAS) Nino3 
and California Nino/Nina indices is —0.65 (—0.32) in the past 30 
years, significant at the 99.99% (90%) confidence level by the two 
tailed t test. If the interannual variations of the California Nino/Nina 
indices simultaneously related to ENSO are linearly regressed out, 
the California Nino/Nina shows the highest interannual variability in 
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Figure 1 | (a) Standard deviation and (b) the first EOF mode of the monthly SST anomalies and (c) the California Nino/Nina indices and (d) their 
seasonally-stratified standard deviations (grey bar) for the period of January 1982 to December 201 1. The enclosed coastal ocean by the dark 
frame ( 1 10°W-120°W, 20 "N-30 "N) in (a) is the region where the averaged SST anomalies are defined as the indices of California Nino/Nina. The dark 
line in (d) denotes the seasonally-stratified standard deviations of California Nino/Nina indices after the simultaneous variations related to ENSO are 
linearly regressed out. The figure was plotted by GrADs software. 
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Figure 2 | (a) Lead-lag correlation coefficients between the JAS California Nino/Nina indices and 3-month-running mean anomalies in SST (shading), 
SLP (contour) and 10-meter-height wind (vector). Negative (positive) numbers in the top of each panel denote the months that the JAS California 
Nino/Nina indices lag (lead), (b) is the same as (a) except that the correlation coefficients are calculated by the JAS California Nino/Nina indices after 
linearly regressing out the simultaneous variations related to ENSO. Correlation coefficients with SST and wind significant at the 95% confidence level 
(—0.4) by the two-tailed t test are shown only. The figure was plotted by GrADs software. 



boreal summer. Since the mixed-layer depth is shallower and the 
isotherms at the surface and in the upper ocean off Baja California 
are more densely packed in boreal summer than winter (Supple- 
mentary Fig. 1), the highest interannual variability in boreal summer 
may be more closely related to the possible coastal Bjerknes feedback. 
The standard deviation and first EOF mode based on the JAS SST 
anomalies (Supplementary Fig. 2) show spatial patterns almost 
identical to those based on the monthly SST anomalies for the whole 
calendar year (Fig. 1 a-b) except that the first EOF mode of the former 
explains a little higher percentage of total variances. Hence, hereafter, 



we focus only on the California Nino/Nina in boreal summer (JAS) 
and investigate the above hypothesis in this report. 

Results 

As shown in Fig. 2a, the anomalous negative sea level pressure (SLP) 
and cyclonic atmospheric circulation appear off the coast of 
California three months before the peak phase of the California 
Nino. The associated alongshore component of the surface wind 
anomalies is poleward and thus reduces the year-round coastal 
upwelling and offshore Ekman transport, resulting in the significant 
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Figure 3 | As in Fig. 2, but for the vertical-zonal section of geopotential height (contour) and air temperature (shading) at 25 °N. Coefficients with air 
temperature significant at the 95% confidence level (—0.4) by the two-tailed ttest are shown only. The figure was plotted by GrADs software. 



positive SST anomalies in the coastal ocean. Over the anomalous 
warm coastal ocean, air parcels in the planetary boundary layer 
(PBL) are thus anomalously heated (Fig. 3a). This causes negative 
SLP anomalies that maintain/enhance the anomalous cross-shore 
pressure gradient in the lower troposphere and the poleward along- 
shore surface wind anomalies, which in turn raises the coastal SSTs. 
The close coupling between the lower troposphere and the coastal 
ocean during the development phase of California Nino can be 
clearly seen in Fig. 4a, where the lead-lag correlation coefficients 
between SST, alongshore surface wind and upwelling indices are 
plotted. It is shown that at the peak time, the correlation coefficient 
between the upwelling (also the alongshore wind) and California 
Nino/Nina indices amounts to —0.7, demonstrating clearly that 



the coastal Bjerknes feedback contributes to about a half of the total 
variances of the California Nino/Nina in boreal summer. The 
California Nino/Nina may be also influenced by the tropical ENSO 
in boreal summer. Hence, significant SST anomalies can be seen in 
the eastern tropical Pacific as well in the correlation figures (Fig. 2a). 
However, we stress that the close regional air-sea coupling off Baja 
California during the developing phase exists clearly even after the 
ENSO influences are linearly removed (Figs. 2-4b). 

The climatological equatorward surface winds are weakened dur- 
ing the developing phase of the California Nino, but possible con- 
tributions from the latent heat flux anomalies to the positive SST 
anomalies averaged in the region of our interest are limited 
(Supplementary Figs. 3-4). This is probably because of the spatial 
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Figure 4 | (a) Lead-lag correlation coefficients between the JAS California Nino/Nina indices (JAS CA) and the 3-month-running mean California 
Nino/Nina (CA, grey filled bar), along-shore surface wind (ASW, dark open bar), upwelling (UMI, blue line) and Nifio3 (red line) indices. ASW is 
positive equatorward. Negative (positive) numbers in the x-axis denote the months that the JAS CA lag (lead), (b) as in (a) except that the correlation 
coefficients are calculated by the JAS CA after linearly regressing out the simultaneous variations related to ENSO (JAS CA-ENSO). Correlation 
coefficients of —0.4 are significant at the 95% confidence level by the two-tailed t test. The figure was plotted by GrADs software. 



inconsistence of the latent heat flux anomalies. In the inshore ocean, 
the anomalously warm SSTs due to oceanic processes may offset or 
even surpass the effect of the reduced wind speed on evaporation 
through changes of the near-surface specific humidity; the latent heat 
flux anomalies there are insignificant until two months before the 
peak, and then become significantly increased from the inshore 
ocean to the overlying atmosphere. In contrast, the latent heat flux 
anomalies contribute significantly to the positive SST anomalies in 
the further offshore ocean particularly at the southern edge of the 
SLP anomalies where the anomalous cyclonic surface winds weaken 
the climatological anticyclonic surface winds and thus reduce the 
latent heat fluxes from the ocean to the atmosphere (Fig. 2, 
Supplementary Figs. 1, 4). In addition, the subtropical ocean is 
known to be covered by low stratus clouds due to trade wind inver- 
sion related to subsidence of the subtropical high 23 . The SST-stratus 
cloud feedback, through which the warmer-than-normal SST 
reduces the amount of stratus clouds, increases the downward short- 
wave radiation and thus enhances the warm SST anomalies, may also 
contribute to the positive SST anomalies during the developing phase 
of the California Nino 24 (Supplementary Fig. 4). However, the surface 



thermal forcing as a whole may play a secondary role compared to 
the oceanic processes on the SST anomalies related to California 
Nino during the developing phase; the highest correlation coefficient 
between the net surface heat flux anomalies averaged over the region 
of our interest and the summer California Nino/Nina indices is —0.4 
at most when the California Nino/Nina indices lag two months for 
the case where the linear ENSO influences are excluded 
(Supplementary Fig. 3b). In contrast, during the decaying phase, 
the accumulated positive SST anomalies off Baja California and 
California are dampened quickly by the surface heat fluxes, in par- 
ticular the latent heat fluxes, to the overlying atmosphere as shown by 
much higher correlation coefficients with the California Nino/Nina 
indices (Supplementary Figs. 3-4). 

The anomalous oceanic advection may also contribute to the 
development of California Nino 3 . The poleward alongshore surface 
wind anomalies off Baja California during the developing phase of 
California Nino may reduce the southward transportation of cold 
water by the California Current, i.e., the equatorward wind-driven 
surface current 25 . The wind anomalies decrease the offshore Ekman 
transport and increase the coastal sea surface heights (SSHs) 
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(Supplementary Fig. 5). The anomalous cross-shore gradient of SSHs 
thus formed is consistent with the poleward surface current anom- 
alies. These processes are comparable to the zonal advective feedback 
in the tropical Pacific during the ENSO events 26 . Hence, in a broad 
sense, the coastal ocean is analogous to the tropical ocean as pointed 
out first by Japanese oceanographer Kozo Yoshida and Chinese 
oceanographer Han- Lee Mao in collaboration about a half century 
ago at the Scripps Institution of Ocenography 27,28 when they found a 
similarity between the equatorial upwelling and the coastal upwelling 
as the ocean responses to the surface winds. Here we have extended 
the similarity to ocean-atmosphere coupled processes. 

Discussion 

Although the generation mechanism of the California Nino/Nina in 
boreal summer is quite similar to that of the tropical ENSO, to 
recognize differences is important, too. The former has much smaller 
temporal and spatial scales; it is confined to the coastal ocean and 
lasts only for 1-2 seasons (Figs. 2-3). This is probably because, 
besides the coastal upwelling, the shallow mixed-layer and sharp 
vertical and horizontal temperature gradients in the upper ocean 
are required for the coastal Bjerknes feedback to operate effectively. 
Therefore, the California Nino/Nina generated by the coastal 
Bjerknes feedback is seasonally locked to boreal summer within 
the coastal ocean in the subtropical transition zone. Also, we note 
that the nature of California Nino/Nina in boreal summer on which 
we have focused in this study differs from that of California Nino/ 
Nina in boreal winter 2 " 5,19 . The latter is generally referred to the 
coastal ocean warming/cooling in the northeast Pacific along almost 
the entire western coast of North America. The broadly-extended 
coastal SST anomalies are attributed mainly to the external forcing of 
the tropical ENSO via both oceanic and atmospheric teleconnections 
and/or the ENSO-independent basin-scale atmospheric circulation 
anomalies in mid-latitudes related to the Aleutian Low. Even if the 
coastal Bjerknes feedback discussed here plays a role, it must be much 
less. Actually, the correlation coefficient between January-March 
California Nino/Nina and alongshore wind (upwelling) indices as 
defined in this study is only about —0.27 ( — 0.16). 

Since the coastal Bjerknes feedback can operate theoretically in the 
subtropical coastal oceans in summer, we expect that similar coastal 
Nino/Nina phenomena may exist along the eastern boundaries of the 
subtropical oceans under subtropical highs in both Northern and 
Southern Hemispheres. In the austral summer of 2010/11, extreme 
coastal warming took place off the western coast of Australia and 
devastated dramatically the marine ecosystem. This anomalous 
warming was called Ningaloo Nino 29 . Although it was attributed 
mainly to the oceanic and atmospheric teleconnections related to 
the co-occurring La Nina in the tropical Pacific, it may have been 
amplified by the coastal Bjerknes feedback 18,2930 . 

The California Nino/Nina lasts more than one season and may 
certainly impact the local marine ecosystem. For example, the 
spawning and meridional migration of Pacific sardine in the 
California Current Ecosystem may be influenced by the related 
SST anomalies 31 , and its abundance may be impacted by the upwel- 
ling-induced anomalies in the biological production 32 . Also, the 
California Nino/Nina may increase/decrease the summer precipita- 
tion in Baja California and surface temperature in the nearby coastal 
continent (Supplementary Fig. 6). Hence, accurate seasonal predic- 
tion of the anomalous events could be beneficial to the industrial 
management and associated economic activities. Since the California 
Nino/Nina in boreal summer is less influenced by ENSO than that in 
boreal winter, its seasonal prediction is challenging. So far, the sea- 
sonal predictability of subtropical SST anomalies is mostly originated 
from their robust relationship with the tropical ENSO 33-35 . For 
instance, the 2010/11 Ningaloo Nino can be successfully predicted 
nine months ahead in a coupled ocean- atmosphere model by virtue 
of the co-occurring La Nina 35 . We believe that efforts to predict 



regional intrinsic climate modes such as California Nino/Nina in 
boreal summer may contribute to enhancing skills of seasonal pre- 
diction for various societal applications as well as to deepening our 
knowledge of regional climate variations in the extratropics. 

Methods 

We use the National Oceanic and Atmospheric Administration (NOAA) Optimum 
Interpolation Sea Surface Temperature (OISST) version 2, Pacific Fisheries 
Environmental Laboratory (PEFL) derived upwelling indices, and National Centers 
for Environmental Prediction (NCEP)/National Center for Atmospheric Research 
(NCAR) reanalysis 1 for the period of January 1982 to December 2011. The California 
Nino/Nina index is denned as the monthly SST anomalies averaged over 110 U W- 
120°W and 20°N-30°N (the enclosed coastal ocean in Fig. la). The Nino3 index is the 
monthly SST anomalies averaged over 90" W-150 C W and 5°S-5°N. The JAS 
California Niho/Niha index with variations related to ENSO linearly regressed out 

(K) is computed by the formula: K — Ki — r 12 — - K 2 . Here, Ki is the JAS California 

T>i 

Nino/Nina index, K 2 the JAS Nino3 index, r 12 the correlation coefficient between 
them, D x the standard deviation of Ki, and D 2 the standard deviation of K 2 . The 
alongshore surface wind index is the alongshore component of the 10-meter-height 
wind anomalies averaged over the same ocean as the California Nino/Nina index after 
projecting the surface wind along the coast of Baja California (321' to the north). The 
PEFL upwelling indices are computed from monthly-mean pressure fields and reflect 
the wind-induced coastal upwelling and off-shore Ekman transport at 15 standard 
locations along the west coast of North America 21 . Among them, three locations are 
within the coastal ocean of our interest (24°N 1 13°W, 27°N 1 15°W and 30°N 1 19°W) 
and thus their monthly anomalies are averaged to represent the upwelling index in 
this study. 

Three-month-running mean anomalies are adopted to represent the monthly 
anomalies in Figs. 2-4 and Supplementary Figs. 3-6, 8-10 to minimize the intra- 
seasonal variations, but similar results can be obtained if monthly data are used 
(Supplementary Fig. 7). The linear correlation and regression analyses are applied in 
this study to demonstrate the existence of coastal Bjerknes feedback for simplicity. 
However, there may be asymmetry between California Nino and Nina. For instance, 
California Nina may have larger anomalous amplitudes and spatial extension than 
Nino and last for a longer time (Supplementary Figs. 8-10). Also, the JAS California 
Niho/Niha indices show a negative skewness of —0.11/ — 0.20 with/without the 
simultaneous variations linearly related to ENSO. The asymmetry may be partly due 
to the co-occurrence of negative phase Pacific Decadal Oscillation with California 
Nifia (Supplementary Fig. 8b) and partly due to non-linear processes; this needs to be 
studied further. We note that the results shown in this study do not change quanti- 
tatively even if the NOAA Extended Reconstructed Sea Surface Temperature (ERSST) 
version 3b, NCEP-National Energy Research Supercomputing Center of the 
Department of Energy (DOE) reanalysis 2 and European Centre for Medium-Range 
Weather Forecast Interim Reanalysis (ERA-Int) are used. 
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